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Abstract  In  this  study ,  we  used  the  Navy  ’  s  Generalized  Digital  Environmental  Model 
(GDEM)  climatological  temperature  and  salinity  data  on  a  0.5°  x  0.5°  grid  to  investigate  the  sea¬ 
sonal  variabilities  of  the  southwest  Philippines  Sea  (0.5°  -  9°N,  123.5°  -  136.5°)  thermohaline 
structure  and  circulation .  The  GDEM  for  the  area  was  built  up  on  historical  ( 1930  —  1997 )  temper¬ 
ature  and  salinity  profiles.  A  three-dimensional  estimate  of  the  absolute  geostrophic  velocity  field  on 
isopycnal  surface  was  obtained  from  the  GDEM  temperature  and  salinity  fields  using  the  P-vector 
method .  The  seasonal  variabilities  of  the  thermohaline  structure  and  currents  (  obtained  from  the  in¬ 
verse  method)  such  as  the  Mindanao  Current,  Mindanao  Undercurrent,  North  Equatorial  Counter 
Current,  New  Guinea  Coastal  Undercurrent,  and  dual-eddies  (cyclonic  Mindanao  Eddy  and  anticy- 
clonic  Halmahera  Eddy)  are  identified. 
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1.  INTRODUCTION 

The  western  equatorial  Pacific,  particularly  the  southernmost  Philippine  Sea,  was  called  “water 
mass  crossroads”  by  Fine  et  al.  (1994)  due  to  the  confluence  there  of  several  water  masses  from 
higher  latitudes  of  both  hemispheres  (Wyrtki,  1961 ;  Fine  et  al. ,  1994) .  Fine  et  al.  (1994)  de¬ 
picted  (  Fig .  1 )  major  currents  in  the  Indonesian  region . 

After  encountering  the  western  boundary  along  the  Philippine  coast,  the  North  Equatorial  Cur¬ 
rent  (NEC)  bifurcates  into  the  northward  flowing  Kuroshio  and  the  southward  flowing  Mindanao  Cur¬ 
rent  (MC)  (Nitani,  1972).  The  MC  and  the  New  Guinea  Coastal  Undercurrent  (NGCUC)  flow 
equatorward,  feeding  the  North  Equatorial  Counter  Current  (NECC)  ,  the  New  Guinea  Coastal  Cur¬ 
rent  (NGCC)  ,  and  the  flow  from  the  Pacific  to  the  Indian  Ocean,  that  is,  the  Indonesian  Through- 
flow  ( IT)  .  The  MC  is  unique  because  it  is  the  only  northern  hemisphere  low- latitude  western  bound¬ 
ary  current  that  has  a  mean  flow  toward  the  equator  (Lukas  et  al.  ,  1996)  ,  and  is  important  because 
it  feeds  to  IT  (Ffield  and  Gordon,  1992),  and  impacts  on  the  stratification  of  the  Indian  Ocean 
(Lukas  et  al. ,  1996;  Godfrey,  1996) . 
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Fig .  1  Map  of  the  western  tropical  Pacific  Ocean  and  Indonesian  Seas  showing  the  major  geographic  names  and  surface 
to  intermediate  depth  currents .  Here  ME  and  HE  represent  Mindanao  and  Halmahera  eddies  (  after 
Fine  et  al.  ,  1994) 

The  NECC  divides  into  three  branches  near  130°E.  The  central  branch  continues  to  flow  east¬ 
ward  as  the  NECC .  The  northern  branch  moves  northwestward  and  forms  a  cyclonic  eddy  near  Mind¬ 
anao  Island,  and  is  called  the  Mindanao  Eddy  (ME).  The  southern  branch  moves  southwestwaid ,  joins 
the  northwest  moving  currents  from  NGCC  and  NGCUC ,  which  are  parts  of  the  South  Equatorial  Current 
(SEC),  and  forms  an  anticyclonic  eddy  called  Halmahera  Eddy  (HE)  near  the  Halmahera  Sea. 


Fig. 2  Geography  and  isobaths  showing  the  bottom  topography  of  the  southwestern  Pacific  near  Mindanao  Island 
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Metzger  and  Hurlburt  (1996)  adopted  a  high  horizontal  resolution  six- layer  model  driven  by 
monthly  climatological  winds  to  analyze  the  coupled  dynamics  among  the  Sulu  Sea ,  the  South  China 
Sea,  and  the  Pacific  Ocean.  Their  model  simulates  the  existence  of  the  dual  eddy  (ME  and  HE) 
structure  in  the  middle  of  the  “water  mass  crossroads”  surrounded  by  complicated  topography  (Fig. 
2),  and  connected  to  major  equatorial  currents.  Seasonal  variabilities  (thermohaline  structure  and 
circulation)  of  the  two  eddies  represent  the  seasonal  variability  in  the  low  latitude  dynamics  and  in¬ 
terbasin  water  mass  exchange . 

Water  mass  characteristics  had  been  used  to  infer  flow  along  isopycnal  surfaces.  For  example, 
between  300  and  1000  m  the  North  Pacific  Intermediate  Water  (NPIW)  (around  ae  =26.8)  and 
Antarctic  Intermediate  Water  (AAIW)  (around  Og  =  27 .2)  reach  this  area  (Reid,  1965;  Tsuchiya, 
1991;  Tally,  1993;  Bingham  and  Lukas ,  1994,  1995;  Fineetal.,  1994).  Kashinoetal.  (1996) 
confirmed  the  existence  of  NPIW  and  AAIW  in  the  far  western  equatorial  Pacific  near  the  entrances 
to  the  Celebes  and  Halmahera  Seas .  They  also  provided  further  evidence  for  the  flow  of  waters  from 
the  southern  hemisphere  across  the  equator ,  turning  eastward  into  NECC . 

To  obtain  current  signal  quantitatively  from  hydrographic  data  in  this  area  is  not  easy .  Godfrey 
(1996)  pointed  out  that  there  are  a  number  of  special  difficulties  such  as  active  internal  tides  and 
the  fact  that  the  geostrophy  can  only  be  used  with  caution  to  estimate  the  current  distribution  in  the 
area.  To  overcome  these  difficulties,  we  used  the  P- vector  method  on  the  isopycnal  surfaces  (Chu, 
1995  ,  2000;  Chu  and  Li,  2000;  Chu  et  al. ,  1998,  2001a, b)  to  invert  the  absolute  velocity  from 
the  Navy  ’  s  public  domain  Global  Digital  Environmental  Model  ( GDEM )  climatological  monthly 
mean  temperature  and  salinity  data  set  with  1/2°  x  1/2°  resolution  (Teague  et  al.  ,  1990)  .  The  P- 
vector  method  is  a  two-step  inverse  method  for  determining  first  the  direction  of  the  current  (  P- vec¬ 
tor)  and  then  the  magnitude  using  the  thermal  wind  relation .  The  two  necessary  conditions  for  valid¬ 
ity  of  inversion  are  easily  incorporated  into  the  P- vector  method  (Chu  et  al.  ,  1998)  . 

We  investigated  the  seasonal  variabilities  of  the  three-dimensional  circulation  and  thermohaline 
structure  of  the  dual  eddies  (ME  and  HE)  .  The  rest  of  the  paper  is  outlined  as  follows:  Section  2 
describes  the  GDEM  T ,  S  data  in  the  isopycnal  coordinate  system.  Section  3  depicts  the  monthly 
mean  velocity  fields  calculated  using  the  P-vector  method .  Section  4  discusses  the  seasonal  variabili¬ 
ty  of  the  volume  transport  streamfunction .  Section  5  presents  the  conclusion . 

2.  ISOPYCNAL  COORDINATE  SYSTEM 


Vertical  discretization 


The  GDEM  T,  S  dataset  showed  that  the  potential  density  a9  varies  from  22 .2  to  27 .725 .  We 
discretized  ae  with  the  increment  Aa  =  0.025  kg/m3  to  obtain  222  a9 -layers  within  which  the  density 
is  vertically  uniform .  To  resolve  well  the  isopycnal  surfaces ,  we  used  cubic  spline  to  interpolate  the 
T,  S  data  into  246  z- levels  with  three  different  increments:  5  m  from  0  to  100  m  depth,  10  m  from 
100  to  1000  m  depths,  20  m  from  1000  to  2500  m  depths,  and  50  m  below  2500  m  depths.  Thus, 
we  built  up  a  high  resolution  z-coordinate  dataset  [  T(  z;-  )  ,  S  (  Zj )  ,  ae  (  Zj  )  ]  .  Notice  that  the  value  of 
27 . 725  is  the  maximum  value  for  <j0  ,  computed  from  the  GDEM  T ,  S  dataset . 

Transformation  of  T,  S  data  from  z -coordinates  to  <re -coordinate  system 

The  transformation  was  done  by  comparing  the  z-coordinate  potential  density  data  Og  (  z, )  with 
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the  discrete  a- values,  a(k)  .  The  two  values  refer  to  different  dimensions.  ae  (  Zj  )  for  the  depth-  *. , 
and  a  (A:).  Thus,  we  use  the  sigma  value  at  the  bottom  of  the  kth  isopycnal  layer,  <sb  (k)  for  the 
comparison , 

ob(k.),  =  -^-[<t(A)  +  o(k  +  1)]  (1) 

The  depth  of  the  bottom  of  the  <j(  A) -layer  is  obtained  by 

if  a(zj)  =  fft(k)  (2) 

and 


W  =  zy  + 


(3) 

where  the  superscript  6  indicates  the  bottom  of  the  Ath  isopycnal  layer.  The  thickness  of  the  /cth 

isopycnal  layer  is  obtained  by 

£  (»)  _  y  -4 

"t  =  2*-l  -  (4) 

After  h[}  is  obtained,  we  may  compute  the  potential  vorticity  on  the  isopycnal  surface,  and  then 
use  the  P-vector  method  to  compute  the  absolute  velocity  (Chu  and  Li,  2000;  Chu  et  al.  ,  2001c)  . 


Bi-monthly  mean  thermohaline  fields  on  <re  =25.0  kg/m3 

In  the  pycnocline  layer  ( 100  to  300  m)  are  salinity  maximum  layers  from  the  North  and  South 
Pacific.  These  water  masses  were  called  North  Pacific  Tropical  Water  (NPTW)  and  South  Pacific 
Tropical  Water  (SPTW)  by  Kashino  et  al.  (1996)  .  It  is  not  our  intention  to  present  a  detailed  cli¬ 
matological  thermohaline  structure  in  this  paper.  We  only  show  some  representative  features  on  ae  = 

25.0  kg/m  (pycnocline)  as  illustration.  This  level  is  generally  considered  to  be  where  the  salinity 
maximum  occurs  (Kashino  et  al.  ,  1996)  . 

Temperature  The  bi-monthly  mean  temperature  on  ae  =  25 . 0  (Fig.  3)  shows  a  cool 

(northern) -warm  (southern)  eddy- like  structure  with  strong  seasonal  variability.  The  cool  (warm) 
eddy  is  enclosed  by  the  18.25  (200 )  isoline.  The  cool  eddy  shows  up  almost  all  year  round  and  is 
located  at  128°  -  133°E  and  6°  -  9°N ,  and  the  warm  eddy  is  evident  only  in  winter  (February)  and 
spring  (April)  .  A  weak  thermal  front  exists  between  the  cool  and  warm  eddies  with  horizontal  tem¬ 
perature  gradient  of  less  than  I'd/ 100  km.  This  front  strengthens  from  December  to  April,  and 
weakens  from  April  to  December. 

Salinity  Salinity  is  generally  used  for  identifying  the  water  masses  in  the  area.  For  ex¬ 

ample,  Kashino  et  al.  (1996)  used  the  35.0  psu  contour  as  the  boundary  of  high  salinity  SPTW  in 
analyzing  data  collected  by  the  R/VKaiyo  during  the  World  Ocean  Circulation  Experiment  (WOCE) 
expedition  I  in  October  1992  and  II  in  February  1994 .  Both  expeditions  covered  the  same  area  as 
ours  (Fig. 4)  .  Horizontal  distribution  on  <j9  =  25 .0,  observed  during  Kaiyo  WOCE  I  and  II,  showed 
the  high-salinity  SPTW  in  the  southern  area.  In  October  1992,  SPTW  reached  2°N  near  Morotai  Is¬ 
land  and  3°N  along  134°E  (Fig. 4a)  .  In  February  1994,  SPTW  reached  farther  north  than  in  Octo¬ 
ber  1992  (Fig. 4b) ,  i.e. ,  north  of  Morotai  Island  and  5°N  along  130°E  (Kashino  et  al.  ,  1996) . 

The  monthly  mean  salinity  on  a9  =  25 . 0  (  Fig .  5  )  shows  feature  similar  to  that  depicted  by 
Kashino  (1996) .  The  high-salinity  SPTW  (35.0  psu)  is  located  in  the  southern  part  of  the  area, 
and  expands  toward  the  northwest  in  Feb.-Jun.  ,  and  retreats  toward  the  southeast  in  August-De- 
cember.  SPTW  has  the  largest  penetration  in  February  -  April.  The  boundary  of  SPTW  is  at  2°  -  3° 
N  near  Morotai  Island  (128°E)  ,  reaches  northernmost  (5°N)  at  131°E,  and  southward  turns  to  3°N 
at  136°E.  In  February  and  April,  a  high  salinity  center  enclosed  by  35.3  psu  occurs  at  131°  -  133° 
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3 .  VELOCITY  FIELD 

The  P-vector  method  for  the  isopycnal  coordinate  system  (Chu  and  Li,  2000;  Chu,  2001)  was 
used  to  determine  the  absolute  velocity. 

Horizontal  distribution 

Most  important  features  from  our  computation  are  the  southward  flowing  MC ,  the  eastward  flow¬ 
ing  NECC,  the  northwestward  flowing  NGCC  and  NGCUC  and  associated  mesoscale  eddies  such  as 
the  cyclonic  Mindanao  Eddy  (ME)  ,  and  the  anticyclonic  Halmahera  Eddy  (HE)  .  Fig. 6  shows  the 
horizontal  velocity  vectors  on  three  different-levels,  pycnocline  (as  =25.0  kg/m3),  intermediate  ( ae 
=  26.5  kg/m3 ) ,  and  deep  ( ae  =27.2  kg/m3 )  levels .  The  last  two  levels  may  represent  NPIW  and 
AAIW  in  this  area. 


124  126  128  130  132  134  E 136  124  126  128  130  132  134  136  E 

Fig. 5  Bi-monthly  mean  salinity  (psu)  fields  on  oa  =25.0  kg/m3  from  GDEM 


Pycnocline  level  ae  =25.0  kg/m3  Three  major  currents  (MC,  NGCUC,  and  NECC) 

and  dual  eddies  ( ME  and  HE )  are  easily  identified .  They  have  strong  seasonal  variabilities .  MC 
strengthens  from  October  to  February  and  weakens  from  April  to  August .  After  leaving  the  south  tip 
of  Mindanao  Island,  MC  flows  southeastward  following  the  continental  slope  (Figs  12  and  2)  and  re¬ 
circulates  northeastward  near  Morotai  Island .  Such  a  recirculation  leads  to  the  formation  of  a  cyclonic 
eddy,  ME.  The  MC  is  very  weak  from  April  to  August,  strengthens  in  October,  and  becomes  strong 
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in  the  winter  with  a  maximum  speed  of  around  0.25  m/s  in  February.  NGCUC  flows  northwestward 
along  northern  New  Guinea  Island  ( Fig .  1 )  .  After  leaving  the  north  tip  of  New  Guinea  Island ,  it 
continues  to  move  northward  and  feeds  an  anticyclonic  eddy ,  HE .  This  eddy  strengthens  in  summer 


Fig. 6  Bi-monthly  mean  velocity  vector  field  east  of  Mindanao  Island  on  the  surface  of  (a)  o9  =  25 .0  kg/m3  ,  (b)  o9  = 
26.5  kg/m3  ,  and  (c)  o9  =27.2  kg/m3 
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Fig. 7  Seasonal  variation  of  north-south  velocities  (unit,  cm/s;  northward  positive)  along  7°15'N  latitude  repre¬ 
senting  the  flow  east  of  Mindanao  Island,  (a)  on  ae  levels,  and  (b)  on  z  levels 

Intermediate  level  (ae  =26.5  kg/ m3  )  As  the  water  depth  increases ,  the  seasonal  variation 
of  the  velocity  reduces .  The  intermediate  level  (ae  =  26.5)  is  close  to  the  level  (ae  =26.8  kg/m3 ) 
where  NPIW  located.  The  MC  is  an  energetic,  nearly  150  km  wide,  southward,  coastal  trapped  jet. 
It  turns  east  after  leaving  Mindanao  Island  and  feeds  into  the  NECC ,  which  separates  into  three 


(June  -  August)  and  weakens  in  winter  (December)  .  In  August,  the  HE  is  located  at  128°  -  133° 
E ,  2°  -  6°N  with  a  maximum  tangential  velocity  of  around  30  cm/s . 
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branches  near  130°E.  The  north  and  south  branches  feed  into  two  eddies,  cyclonic  eddy  (ME)  and 
anticyclonic  eddy  (HE)  .  The  central  branch  continues  as  the  NECC .  In  the  winter  (  December-Feb- 
ruary),  the  dual  eddies  are  north-south  oriented  with  ME  (maximum  tangential  velocity  around  15 
cm/s  located  at  128°  -  133°E,  6°  -  9°N  and  HE  (maximum  tangential  velocity  around  25  cm/s  lo¬ 
cated  at  128°  -  133°E,  0.5°  -  4°N.  In  April,  ME  (HE)  expands  toward  the  southwest  (northeast)  . 

In  June ,  the  dual  eddies  are  east- west  oriented  with  ME  in  the  west  and  HE  in  the  east .  The  size  of 
both  eddies  are  smaller  in  the  summer  than  in  the  winter . 

Deep  level  (oe  =  27 .2  kg/m3 )  At  the  deep  level  (<J0  =  27  .2  kg/m3 )  where  AAIW  locat¬ 

ed,  the  flow  pattern  is  quite  different  from  that  on  the  pycnocline  and  intermediate  levels.  An  evi¬ 
dently  westward  flowing  current  is  identified  between  2°  -  6°N  which  is  in  opposite  direction  to  that 
of  NECC  on  the  pycnocline  and  intermediate  level .  The  flow  pattern  clearly  shows  NGCUC  entering 
the  eastern  part  of  the  southern  boundary  (134°  -  136°E)  and  flowing  northwestward  with  a  maxi¬ 
mum  velocity  larger  than  20  cm/s  in  the  winter  (December-February)  . 

In  the  winter  (December  -  February),  NGCUC  bifurcates  at  133° E,  4°N  into  two  branches 
(northward  and  westward)  .  The  northward  branch  recirculates  at  6°N  and  forms  an  anticyclonic  me- 
soscale  eddy  (around  200  km  in  diameter)  .  The  westward  branch  flows  zonally  (under  NECC)  to  a 
longitude  near  129°E ,  and  splits  into  three  branches  moving  southward ,  westward ,  and  northward . 
The  westward  branch  transports  AAIW  into  the  Celebes  Sea.  The  northward  branch  turns  to  west  at 
8°N,  recirculates  along  the  west  coast  of  Mindanao  Island,  and  forms  a  cyclonic  eddy  (i.e.  ,  ME)  . 
In  the  summer  (June  -  August)  ,  a  200  km  diameter  cyclonic  eddy  with  maximum  swirl  velocity  of 
15  cm/s  is  identified  southeast  of  the  southern  tip  of  Mindanao  Island . 

Major  currents 

MC  (along  7°15’N)  As  we  mentioned,  the  MC  is  a  major  current  at  the  water  mass 

crossroads  and  feeds  to  NECC  and  IT.  The  seasonal  variability  of  MC  affects  the  seasonal  variability 
of  NECC  and  IT,  which,  in  turn,  impact  on  the  seasonal  variability  of  ME  (due  to  the  NECC  vari¬ 
ability)  and  the  Indian  Ocean  stratification  (due  to  the  IT  variability)  . 

We  use  the  north-south  geostrophic  velocities  on  isopycnal  surfaces  along  7°15’N  to  represent 
the  seasonal  variability  of  MC  (Fig. 7a).  Positive  (negative)  values  show  the  poleward  (equator- 
waid)  flow.  The  negative  values  (  equatorward  flow)  are  shaded  by  grey  color.  The  top-to-bottom 
black  vertical  bars  represent  the  islands .  The  small  black  bars  at  the  top  or  the  bottom  represent  that 
ae  does  not  reach  the  values  there.  The  equatorward  flowing  MC  has  two  jet-cores,  located  in  the 
upper  layer  above  the  pycnocline  (from  the  surface  to  a9  =  25.0)  ,  and  deep  layer  (a9  >  26.5)  ,  re¬ 
spectively  .  The  upper  layer  jet-core  has  strong  seasonal  variability .  Evident  equatorward  MC  occurs 
from  December  to  March  at  a  maximum  speed  of  10  cm/s  in  January  and  February  and  5  cm/s  in 
March  near  the  east  coast  of  Mindanao  Island  at  a9  =23.25.  MC  disappears  from  April  to  Novem¬ 
ber.  The  deep  layer  jet-core  exists  with  weak  seasonal  variability  all  year  round .  Its  strength  reaches 
the  maximum  (  >  25  cm/s)  in  the  winter  (February)  and  the  minimum  in  the  summer  (June  -  Au- 
gust)  . 

We  convert  the  ^-component  on  o9  to  that  on  z  levels  (Fig. 7b)  .  The  ^-component  on  z-level 
also  shows  the  upper  layer  (around  50  m)  jet-core  and  the  deep  layer  jet-core  (300  -  600  m)  ,  re¬ 
spectively.  A  maximum  speed  of  100  cm/s  in  MC  was  measured  along  7  N  (Lukas  et  al.  ,  1991) 
with  the  use  of  Acoustic  Doppler  Current  Profilers  (ADCP)  .  The  calculated  absolute  velocities  from 
GDEM  underestimate  the  actual  velocities . 
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Fig. 8  Seasonal  variation  of  north- south  velocities  (unit:  cm/s;  northward  positive)  along  1°15'N  latitude  representing 
the  flow  east  of  Halmahera  Island:  (a)  ono,  levels,  and  (b)  on  z  levels 

Fig .  7b  shows  the  existence  of  three  northward  velocity  cores .  Among  them  two  cores  are  close 
to  the  coast  (west  of  130°E)  and  one  core  (centered  at  300  m  depth)  is  away  from  the  coast  (130° 
-  134°E)  .  The  two  coastal  cores  represent  the  Mindanao  Undercurrent  (MUC)  described  by  Hu  and 
Cui  (1989)  ,  and  the  third  core  represents  the  eastern  part  of  the  ME. 
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The  MUC  has  evident  seasonal  variation.  From  April  to  September,  the  upper  layer  MUC,  oc¬ 
curring  from  50  to  300  m  (centered  at  200  m)  ,  is  close  to  the  coast  and  has  a  width  of  around  200 
km  and  a  maximum  velocity  of  10  cm/s .  It  retreats  toward  the  coast  from  October  to  December  with 
a  width  less  than  100  km,  and  expands  toward  the  east  to  a  longitude  near  132°E  from  January  to 
March  with  a  maximum  velocity  of  10  cm/s  at  around  150  m  deep.  The  lower  layer  MUC,  occurring 
below  400  m,  is  around  100  km  away  from  the  coast  and  has  a  weak  seasonal  variability  in  its  loca¬ 
tion  and  a  strong  seasonal  variability  in  its  strength.  It  reaches  a  maximum  strength  (20  cm/s)  in 
February  and  a  minimum  strength  (5  cm/s)  from  March  to  October. 

ME  (Along  7°  15’ N)  The  ME  is  a  sub-surface  system,  represented  by  the  alternate 

southward  flow  west  of  and  northward  flow  east  of  130°E  (Fig. 7b)  .  In  November,  the  ME  extends 
from  200  to  500  m  with  a  jet  core  at  300  m  depth .  From  November  to  March ,  its  position  keeps 
quite  steady ,  but  its  strength  varies  drastically .  The  ME  velocity  enhances  through  February  and 
weakens  in  March.  The  northward  branch  velocity  increases  to  20  cm/s  in  December  and  January, 
and  the  southward  branch  velocity  increases  to  25  cm/s  in  February .  From  April  to  October,  the  ME 
shifts  southwestward  to  128°  —  133°E  and  south  of  6°N  (Fig. 6) . 

NGCUC  (Along  1°15’N)  Similar  to  the  MC,  the  NGCUC  is  a  major  current  flowing 

from  the  Southern  Hemisphere  to  the  water  mass  crossroads  and  feeds  NECC  and  IT .  The  seasonal 
variability  of  NGCUC  affects  the  seasonal  variability  of  NECC  and  IT,  and,  in  turn,  impacts  on  the 
seasonal  variability  of  ME  (due  to  NECC  variability)  and  the  Indian  Ocean  stratification  (due  to  IT 
variability)  . 

Since  the  geostrophic  balance  is  not  generally  good ,  we  use  the  computed  v  velocities  on  isopy- 
cnal  surfaces  along  1°15’N  to  calculate  the  seasonal  variability  of  NGCUC  approximately  (Fig. 8a)  . 
Positive  (negative)  values  show  the  poleward  (equatonvard)  flow.  The  negative  values  (equatorward 
flow)  are  shaded  by  grey  color.  Same  as  in  Fig. 7a,  the  top-to-bottom  black  vertical  bars  represent 
the  islands.  The  vertical  bars  between  127°15'  and  128°45'E  represent  Halmahera  Island  (Fig. 2)  . 
The  small  black  bars  at  the  top  or  the  bottom  represent  that  <Ja  does  not  reach  the  value  there .  The 
NGCUC  has  two  branches ;  a  400  km  wide  coastal  branch  along  the  east  coast  of  Halmahera  Island 
and  an  around  300  km  wide  open  water  branch  east  of  133°E. 

An  interesting  feature  is  the  out-of-phase  variation  between  the  coastal  and  open  water  branch¬ 
es.  For  example,  the  coastal  (open  water)  branch  strengthens  (weakens)  from  September  to  March, 
and  the  coastal  (open  water)  branch  weakens  (strengthens)  in  spring  from  March  to  September. 
From  January  to  April ,  the  coastal  branch  is  strong  with  two  cores ,  one  in  the  upper  layer  ( above 
the  pycnocline ,  i .  e . ,  ae  <  25 . 0 )  and  the  other  in  the  intermediate  layer  (  25 . 0  <  a9  <  26 . 5  )  .  The 
maximum  northward  velocity  is  around  15  cm/s.  The  open  water  branch  is  located  in  the  deep  layer 
(ae  >26.5)  and  has  maximum  northward  velocity  of  15  cm/s.  From  April  to  October,  the  coastal 
(open  water)  branch  weakens  (strengthens)  at  a  maximum  northward  velocity  around  5  cm/s  (40 
cm/s  in  August)  .  Besides,  the  open  water  branch  extends  from  near  surface  to  deep  layer.  Such  an 
out-of-phase  variation  is  also  easily  seen  in  the  cross-section  of  the  component  on  z  level  (Fig. 8b)  . 

NECC  (Along  130°E)  The  NECC,  originated  from  the  confluence  of  southward  flowing 

MC  and  northward  flowing  NGCUC,  is  an  eastward  equatorial  flow  centered  at  around  5°N  (Lukas  et 
al. ,  1991)  .  We  use  the  east-west  geostrophic  velocities  on  isopycnal  surfaces  along  130°E  longitude 
to  represent  the  seasonal  variability  of  NECC  (Fig. 9a)  .  Positive  (negative)  values  show  the  east¬ 
ward  (westward)  flow.  The  negative  values  (westward  flow)  are  shaded  in  grey. 

The  NECC  extends  from  the  surface  to  the  level  of  <j„  =  27  all  year  round .  It  occurs  between  3° 
and  8°N  from  June  to  December  and  south  of  8°N  with  two  branches  separated  by  an  upper  layer 
westward  current  from  January  to  May ,  with  two  jet-cores ,  located  in  the  pycnocline  layer  (  centered 
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around  ae  =25.0),  and  in  the  intermediate  layer  (27.0  >  a0  >  26.0) ,  respectively .  The  pycnocline 
layer  jet-core  has  strong  seasonal  variability  and  maximum  speed  of  40  cm/s  (10  cm/s)  in  August 
(February)  .  The  intermediate  layer  jet-core  is  usually  stronger  than  the  pycnocline  layer  core  with 
less  seasonal  variability .  The  maximum  speed  in  the  intermediate  layer  jet-core  is  45  cm/s  in  Febru¬ 
ary  and  20  cm/s  in  May. 


Fig. 9  Seasonal  variation  of  east- west  velocities  (unit,  cm/s;  eastward  positive)  along  130°E  longitude  repre¬ 
senting  NECC.  (a)  onu,  levels,  and  (b)  on  z  levels 
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HE  (Along  130°E)  We  convert  the  v -component  on  ae  to  that  on  z  levels  (Fig. 9b)  . 

The  NECC  separates  near  130°E  into  three  branches,  central  branch  (continuation  of  NECC)  ,  north 
branch  feeding  to  ME ,  and  south  branch  feeding  to  HE .  The  HE  is  seen  by  the  alternate  eastward 
NECC  and  westward  flow  south  of  NECC .  The  strength  of  HE  is  represented  by  the  velocity  of  the 
westward  flow .  The  HE  strengthens  in  summer  ( June  -  August )  and  weakens  in  winter  ( Decem¬ 
ber)  .  The  maximum  westward  speed  in  the  southern  flank  of  the  HE  reaches  20  cm/s  in  June  and 
July  and  10  cm/s  in  February  and  March. 
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Fig.  10  Vertically  integrated  velocity  vectors 


4.  VOLUME  TRANSPORT 


Vertically  integrated  velocity 

Monthly  mean  vertically  integrated  (for  the  whole  water  column)  velocity  is  computed  to  repre¬ 
sent  the  volume  transport  (Fig.  10)  .  Three  major  currents  (MC,  NGCUC,  and  NECC)  and  dual  ed¬ 
dies  (ME  and  HE)  are  also  identified  in  the  integrated  flow  field  whose  seasonal  variation  is  evi¬ 
dent  .  The  integrated  MC  strengthens  from  October  to  February  and  weakens  from  April  to  August . 
After  leaving  the  south  tip  of  Mindanao  Island ,  MC  flows  southeastward  following  the  continental 
slope  (Figs  16,  2)  and  recirculates  northeastward  near  Morotai  Island  and  then  turns  towards  the 
east .  The  recirculation  leads  to  the  formation  of  ME  (  cyclonic  eddy )  ,  and  the  eastward  flow  feeds 
into  NECC.  NGCUC  flows  northwestward  into  the  area  between  134°  -  136°E  at  0.5°N  latitude  and 
then  joins  the  NECC  at  133°E,  3°N.  The  ME-HE  dipole  pattern  is  evident  from  June  to  August . 

MC  volume  transport 

We  computed  monthly  total  and  layered  (between  two  ae  levels)  volume  transports  across  7°15' 
N  (northward  positive)  between  126°45'  (Mindanao  coast)  to  lSO'M^E  longitude  (Fig.  11)  to  rep¬ 
resent  the  seasonal  variability  of  MC  volume  transport.  The  dashed,  dotted,  and  solid  curves  indi- 
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Fig.  12  Monthly  variation  of  total  and  layered  (between  two  ae  levels)  zonal  volume  transports  across  130°E  longitude  between 
(745'  and  6°15'N,  representing  the  NECC  transport  (Sv).  (a)  total,  (b)  o9  less  than  24.0,  (c)  <r9  from 
24.0  to  25.0,  (d)  a9  from  25.0  to  26.0,  (e)  a9  from  26.0  to  27.0,  and  (f)  a9  more  than  27.0.  Hie 
dashed,  dotted,  and  solid  curves  indicate  the  eastward,  westward,  and  net  transports 

Volume  transport  to  the  Celebes  Sea 

We  commuted  monthly  total  and  layered  (between  two  ae  levels)  volume  transports  across  125° 
15'E  (eastward  positive)  between  2°15'  to  5°45'N  latitude  (Fig.  13)  to  represent  the  seasonal  vari¬ 
ability  of  volume  transport  to  the  Celebes  Sea.  The  dashed,  dotted,  and  solid  curves  indicate  east¬ 
ward  ,  westward ,  and  net  transports .  The  eastward  flow  is  much  weaker  than  the  westward  flow  for 
the  total  transport  (denoting  water  entering  the  Celebes  Sea)  in  January  (  -  16.7  Sv)  .  The  strength 
of  the  volume  transport  of  MC  into  the  Celebes  Sea  weakens  from  winter  to  spring,  and  in  summer 
(August),  the  net  transport  (4.6  Sv)  is  eastward,  indicating  outflow  from  the  Celebes  Sea  to  the 
Pacific  Ocean.  The  most  layered  transports  are  negative  (westward)  all  year  round  except  the  mid¬ 
layer  (  26 . 0  ^  <je  sg  25 . 025  )  ,  where  the  eastward  flow  is  dominant . 

Volume  transport  to  the  Molucca  Sea 

We  computed  monthly  total  and  layered  (between  two  ae  levels)  volume  transports  across  1°45'N 
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Month 

Fig.  13  Monthly  variation  of  total  and  layered  (between  two  <J6  levels)  zonal  volume  transports  across  125°15'E  (eastward 
positive)  longitude  between  2°15'  and  5°45'N,  representing  the  volume  transport  (Sv)  between  the  Celebes  Sea  and  Pacific 
Ocean,  (a)  total,  (b)  o9  less  than  24.0,  (c)  o9  from  24.0  to  25.0,  (d)  o6  from  25.0  to  26.0,  (e)  <J9  from  26.0  to  27 .0, 
and  (f)  o9  more  than  27.0.  Hie  dashed,  dotted,  and  solid  curves  indicate  the  eastward,  westward,  and  net  transports 
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Fig .  14  Monthly  variation  of  total  and  layered  ( between  two  ae  levels )  latitudinal  volume  transports  ( northward  positive ) 
across  1°45'N  latitude  between  125°45'  and  127°25'E,  representing  the  seasonal  varia  bility  of  volume  transport  to  the  Moluc¬ 
ca  Sea.  (a)  total,  (b)  <j9  less  than  24.0,  (c)  <j9  from  24.0  to  25.0,  (d)  <J9  from  25.0  to  26.0,  (e)  <J9  from  26.0  to  27 .0, 
and  (f)  <J9  more  than  27.0.  The  dashed,  dotted,  and  solid  curves  indicate  the  northward,  southward,  and  net  transports 
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(northward  positive)  between  125045'  to  127°25’E  longitude  (Fig.  14)  to  represent  the  seasonal 
variability  of  volume  transport  to  the  Molucca  Sea .  The  dashed ,  dotted .  and  solid  curves  indicate 
northward,  southward,  and  net  transports.  The  southward  flow  (into  the  Molucca  Sea)  is  a  little 
weaker  than  the  northward  flow  (out  of  the  Molucca  Sea)  for  the  total  transport  in  January  (1.5 
Sv)  .  The  net  northward  volume  transport  decreases  with  time .  From  February  to  May ,  the  net  vol¬ 
ume  transport  is  southward  with  the  minimum  value  of  -  1 . 8  Sv  in  May .  After  May ,  the  net  volume 
transport  becomes  northward  and  enhances  from  2.8  Sv  in  June  to  7.9  Sv  (maximum  value)  in  Au¬ 
gust.  The  southward  volume  transport  (into  the  Molluca  Sea)  is  evident  from  January  to  June  in  the 
mid-layer  (  26 . 0  ^  ae  ^  25 . 025  )  . 


Month 


Fig.  15  Monthly  variation  of  total  and  layered  (between  two  ae  levels)  latitudinal  volume  transports  (northward  positive) 
across  1°45'N  latitude  between  129°30'  and  136°30'E,  representing  the  seasonal  variability  of  vol¬ 
ume  transport  to  the  Halmahera  Sea.  (a)  total,  (b)  o9  less  than  24.0,  (c)  o9  from  24.0  to  25.0, 

(d)  ae  from  25.0  to  26.0,  (e)  o9  from  26.0  to  27.0,  and  (f)  o9  more  than  27.0.  Hie  dashed, 
dotted,  and  solid  curves  indicate  the  northward,  southward,  and  net  transports 

Volume  transport  to  the  Halmahera  Sea 

We  computed  monthly  total  and  layered  (between  two  <J9  levels)  volume  transports  across  1°45' 
N  (northward  positive)  between  129°30'  to  136°30'E  longitude  (Fig.  15)  to  represent  the  seasonal 
variability  of  volume  transport  to  the  Halmahera  Sea .  The  dashed ,  dotted ,  and  solid  curves  indicate 
northward,  southward,  and  net  transports.  The  southward  flow  (into  the  Halmahera  Sea)  is  stronger 
than  the  northward  flow  (out  of  the  Halmahera  Sea)  for  the  total  transport  in  most  of  the  year  except 
in  June  and  July.  The  net  southward  volume  transport  (into  the  Halmahera  Sea)  reaches  the  maxi¬ 
mum  strength  (  -  19 .6  Sv)  in  January.  The  net  southward  volume  transport  weakens  from  January  to 
March  (  -  4.5  Sv)  and  then  strengthens  from  March  to  May  (  -  16.4  Sv)  with  time.  The  net  vol¬ 
ume  transport  is  positive  (out  of  the  Halmahera  Sea)  in  June  (5.6  Sv)  and  July  (11.9  Sv)  .  The 
net  volume  transport  is  southward  from  August  (  —  3.5  Sv)  to  December  (  —  18.2  Sv) .  The  south- 
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ward  volume  transport  (into  the  Halmahera  Sea)  is  dominant  in  the  deep  layers  (27.0^a9>26. 

025)  and  ( 27 . 725  Ss  a9  5=  27.025). 

5.  CONCLUSIONS 

(1)  The  GDEM  data  on  the  pycnocline  layer  (ae  =  25.0)  shows  the  cool  ( northern) -warm 
(southern)  dipole  eddy- like  structure  and  the  high-salinity  SPTW  (^35.0  psu)  co- located  with  the 
warm  eddy .  The  SPTW  expands  toward  the  northwest  in  February  -  June  with  largest  penetration  in 
February  -  April ,  and  retreats  toward  the  southeast  in  August  -  December .  The  cool  eddy  occurs  at 
128°  -  133°E  and  6°  -  9°N  almost  all  year  round,  and  the  warm  eddy  is  evident  from  February  to 
April  when  the  SPTW  has  the  largest  northward  penetration .  A  weak  thermal  front  exists  between  the 
cool  and  warm  eddies.  This  front  strengthens  from  December  to  April,  and  weakens  from  April  to 
December . 

(2)  Four  major  currents  (Mindanao  Current,  Mindanao  Undercurrent,  New  Guinea  Coastal 
Undercurrent ,  North  Equatorial  Counter  Current )  and  two  major  eddies  (  Mindanao  Eddy ,  and 
Halmahera  Eddy)  and  their  seasonal  variations  are  identified  from  the  computed  velocity  field.  The 
equatorward  flowing  Mindanao  Current  has  two  jet-cores  located  in  the  upper  layer  above  the  pycno¬ 
cline  (from  the  surface  to  ae  =  25.0),  and  deep  layer  (a9  >  26.5) .  The  upper  layer  jet-core  has 
strong  seasonal  variability  and  maximum  speed  of  10  cm/s  in  January  and  February  and  5  cm/s  in 
March  near  the  east  coast  of  Mindanao  Island  at  ae  =23.25.  It  disappears  from  April  to  November . 
The  deep  layer  jet-core  exists  all  year  round  with  a  weak  seasonal  variability .  Its  strength  reaches  the 
maximum  (  >  25  cm/s)  in  February  and  the  minimum  in  June- August.  The  inverse  estimate  of  the 
maximum  southward  component  in  the  Mindanao  Current  ( 25  cm/s )  is  much  smaller  compared  to 
100  cm/s  from  the  ADCP  measurement  (7°N)  in  August  1988. 

(3)  The  northward  flowing  Mindanao  Undercurrent  with  two  jet-cores  is  identified  near  the 
Mindanao  coast  from  the  computed  ^-component.  From  April  to  September,  the  upper  layer  core, 
occurring  from  50  to  300  m  (centered  at  200  m)  ,  is  close  to  the  coast  and  has  a  width  of  around  200 
km  and  a  maximum  velocity  of  10  cm/s .  It  retreats  toward  the  coast  from  October  to  December  with 
a  width  of  less  than  100  km,  and  expands  toward  the  east  to  a  longitude  near  132°E  from  January  to 
March  with  a  maximum  velocity  (around  150  m  depth)  .  The  lower  layer  core,  occurring  below  400 
m,  is  around  100  km  away  from  the  coast  and  has  weak  seasonal  variability  in  its  location  and  strong 
seasonal  variability  in  its  strength.  It  reaches  maximum  strength  (20  cm/s)  in  February  and  mini¬ 
mum  strength  (5  cm/s)  from  March  to  October. 

(4)  Both  the  Mindanao  Eddy  (cyclonic)  and  Halmahera  Eddy  (anticyclonic)  are  identified.  In 
November,  the  Mindanao  Eddy  extends  from  200  to  500  m  with  a  jet  core  at  300  m  depth .  From  No¬ 
vember  to  March ,  its  position  keeps  steady ,  but  its  strength  varies  drastically .  The  northward  branch 
velocity  increases  to  20  cm/ s  in  December  and  January ,  and  the  southward  branch  velocity  increases 
to  25  cm/s  in  February.  From  April  to  October,  the  ME  shifts  southwestward  to  128°  -  133°E  and 
south  of  6°N .  The  Halmahera  Eddy  strengthens  in  Jun .  -  Aug .  and  weakens  in  December .  The  maximum 
westward  speed  in  the  southern  flank  of  HE  reaches  20  cm/s  in  June  and  July  and  10  cm/s  in  February 
and  March. 

(5 )  The  monthly  total  and  layered  ( between  two  as  evels )  volume  transport  were  computed 
along  7°15'N  (northward  positive)  from  the  Mindanao  coast  to  130°45'E  longitude  to  represent  the 
seasonal  transport  variability  of  the  Mindanao  Current,  and  along  130°E  (eastward  positive)  between 
0°45'  and  8°15'N  to  represent  the  seasonal  transport  variability  of  the  North  Equatorial  Counter  Cur- 


No. 3 


CHU  et  al. :  DETERMINATION  OF  CURRENT  SYSTEM  ON  ISOPYCNAL  SURFACE 


211 


rent .  The  total  transport  of  the  Mindanao  Current  is  equatorward  all  year  round  with  a  maximum  val¬ 
ue  of  40 . 2  Sv  in  January  and  a  minimum  value  of  5 . 3  Sv  in  June  and  July  with  the  annual  mean 
transport  of  23 . 4  Sv ,  which  agrees  with  many  existing  estimates .  The  total  transport  of  the  North 
Equatorial  Countercurrent  is  eastward  all  year  round  with  a  maximum  value  of  45 . 7  Sv  in  January 
and  a  minimum  value  of  20 . 0  Sv  in  April  with  the  annual  mean  transport  of  34 . 4  Sv . 

(6)  The  monthly  total  and  layered  volume  transports  were  computed  along  125°15'E  (eastward 
positive)  between  2°15'  to  5°45'N  latitude  to  represent  the  seasonal  variability  of  volume  transport  to 
the  Celebes  Sea .  The  net  volume  transport  is  into  the  Celebes  Sea  in  most  of  the  year ,  and  has  max¬ 
imum  strength  in  January  (  -  16.7  Sv)  .  It  weakens  from  winter  to  spring,  and  in  summer  (Au¬ 
gust)  ,  the  net  transport  (4.6  Sv)  is  out  of  the  Celebes  Sea  into  the  Pacific  Ocean.  Most  layered 
transports  are  negative  (  westward  )  all  year  round  except  the  mid-layer  (  26 . 0  >  <r9  3:  25 . 025  )  , 
where  the  eastward  (out  of  the  Celebes  Sea)  flow  is  dominant. 

(7)  The  monthly  total  and  layered  volume  transports  were  computed  along  1°45'N  (northward 
positive)  between  125°45'  to  127°25'E  longitude  to  represent  the  seasonal  variability  of  volume 
transport  to  the  Molucca  Sea.  The  net  volume  transport  is  into  the  Molucca  Sea  from  February  to 
May  with  maximum  strength  in  May  (  -  1 .8  Sv)  and  out  of  the  Molucca  Sea  the  rest  of  the  year  with 
maximum  strength  in  August  (7.9  Sv)  .  The  volume  transport  into  the  Molucca  Sea  is  evident  from 
January  to  June  in  the  mid-layer  (26 .0><j3  0^25 .025)  . 

(8)  The  monthly  total  and  layered  volume  transports  were  computed  along  1°45'  (northward 
positive)  between  129°  to  13°30'E  longitude  to  represent  the  seasonal  variability  of  volume  transport 
to  the  Halmahera  Sea .  The  net  volume  transport  is  into  the  Halmahera  Sea  in  most  of  the  year  except 
in  June  and  July.  The  net  volume  transport  into  the  Halmahera  Sea  reaches  maximum  strength  (  - 
19.6  Sv)  in  January.  It  weakens  from  January  to  March  (  -  4.5  Sv)  and  then  strengthens  with  time 
from  March  to  May  (  -  16  .4  Sv)  .  The  net  volume  transport  is  out  of  the  Halmahera  Sea  in  June  (5 . 
6  Sv)  and  July  (11.9  Sv)  .  The  net  volume  transport  is  into  the  Halmahera  Sea  again  from  August 
(-3.5  Sv)  to  December  (  -  18 .2  Sv)  .  The  volume  transport  into  the  Halmahera  Sea  is  dominant 
in  the  deep  layers  (27 .0>a9 ^26 .025)  and  ( 27 . 725 ^ a9 ^ 27 . 025 ) . 
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